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Abstract 

Bounds  on  variables  are  often  implemented  as  part  of  a  quality  control  program  to 
ensure  a  sufficient  pedigree  of  a  product  component  and  these  bounds  may  significantly 
affect  the  product’s  design  through  constraints  such  as  cost,  manufacturability  and 
reliability.  Thus,  it  is  useful  to  determine  the  sensitivity  of  the  product  reliability  to  the 
imposed  bounds.  In  this  work,  a  method  to  compute  the  partial  derivatives  of  the 
probability-of-failure  and  the  response  moments,  such  as  mean  and  the  standard 
deviation,  with  respect  to  the  bounds  of  truncated  distributions  are  derived  for  rectangular 
truncation.  The  sensitivities  with  respect  to  the  bounds  are  computed  using  a 
supplemental  “flux”  integral  that  can  be  combined  with  the  probability-of-failure  or 
response  moment  information.  The  formulation  is  exact  in  the  sense  that  the  accuracy 
depends  only  upon  the  numerical  algorithms  employed.  The  flux  integral  is  formulated  as 
a  special  case  of  the  probability  integral  for  which  the  sensitivities  are  being  computed. 

As  a  result,  the  methodology  can  be  implemented  with  any  probabilistic  method,  such  as 
sampling,  first  order  reliability  method,  conditional  expectation,  etc.  Moreover,  the 
maximum  and  minimum  values  of  the  sensitivities  can  be  obtained  without  any  additional 
computational  cost.  The  methodology  is  quite  general  and  can  be  applied  to  both 
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component  and  system  reliability.  Several  numerical  examples  are  presented  to 
demonstrate  the  advantages  of  the  proposed  method.  In  comparison,  the  examples  using 
Monte  Carlo  sampling  demonstrated  that  the  flux-based  methodology  achieved  the  same 
accuracy  as  a  standard  finite  difference  approach  using  approximately  4  orders  of 
magnitude  fewer  samples.  This  is  largely  due  to  the  fact  that  this  method  does  not  rely 
upon  subtraction  of  two  near-equal  numbers. 


Notation 

a  lower  boundary  of  truncated  distribution 
b  upper  boundary  of  truncated  distribution 
E  expected  value  operator 

fx  (a)  marginal  PDF  of  truncated  distribution  evaluated  at  lower  boundary  a 
fx  (b)  marginal  PDF  of  truncated  distribution  evaluated  at  upper  boundary  b 
/X(x)  joint  probability  density  function  (JPDF) 
g  limit  state 

/(x)  indicator  function  ( 7(x)  =  1  in  failure  region,  and  0  otherwise) 

JPDF  joint  probability  density  function  represented  by  /X(x) 

N  number  of  random  variables 

Pf  probability-of-failure 

Pj  flux  integral  of  probability-of-failure  at  bound  0 
X  vector  of  random  variables,  X  =  (X1,...,XJV) 

Y  vector  of  random  variables  X  excluding  Xn  Y  =  (Xp.  ..,XM,Xi+1,.  ..,XN) 
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Z  response  function  of  random  variables 

Tt  flux  of  JPDF  in  failure  domain  across  the  truncation  boundary 

rz  flux  of  response  times  the  JPDF  across  the  truncation  boundary 

Fz 2  flux  of  response  squared  times  the  JPDF  across  the  truncation  boundary 

Ka  kernel  function  with  respect  to  lower  bound 

Kb  kernel  function  with  respect  to  upper  bound 

)lz  mean  of  response  function  Z 

oz  standard  deviation  of  response  function  Z 

6  lower  or  upper  bound  of  PDF 

1.  Introduction 

In  any  design  problem  it  is  important  to  define  approved  allowables  for  material 
properties,  joint  properties,  fracture  properties,  dimensioning  tolerances,  etc.  Oftentimes 
allowables  arise  from  and  are  implemented  within  a  quality  assurance  (QA)  program 
which  may  establish  bounds  on  allowables  and  ensure  that  values  outside  of  the 
allowables  will  not  be  accepted  in  order  to  guarantee  the  quality  of  a  product.  For 
example,  cases  in  point  include  “A”  (99%  probability,  95%  confidence)  and  “B”  (90% 
probability,  95%  confidence)  basis  allowables.  That  is,  for  an  “A”  basis  allowable,  there 
is  a  95  percent  probability  that  99  percent  of  the  samples  will  be  less  than  or  greater  than 
(depending  on  the  application)  the  numerical  allowable.  For  a  “B”  basis  allowable,  there 
is  a  95  percent  probability  that  90  percent  of  the  samples  will  be  less  than  or  greater  than 
the  numerical  allowable. 
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Another  common  occurrence  arises  with  respect  to  geometric  random  variables 
obtained  by  machining.  Tolerance  limits  (plus  and  minus)  are  placed  upon  dimensions 
that  are  the  result  of  machining  operations,  e.g.,  holes,  fillet  radii,  and  seal  diameters. 
Correct  dimensioning  in  the  part  is  necessary  for  subsequent  assembly  and  to  minimize 
stress  concentrations.  The  QA  operations  ensure  that  dimensions  outside  of  the  tolerances 
will  not  be  present  in  the  finished  product. 

For  these  cases,  it  is  easy  to  see  that  the  implementation  of  a  QA  program  will 
often  result  in  truncated  distributions  that  affect  the  overall  quantity  or  cost  of  a  product. 
If  material  properties  are  outside  of  the  allowables  or  dimensions  outside  of  tolerances 
are  removed  before  production,  then  effectively,  the  quantity  can  be  modeled  as  a 
truncated  distribution.  Typical  cases  include  truncated  normal,  Weibull  and  uniform 
distributions. 

Truncated  distributions  can  also  arise  naturally  when  modeling  nonparametric 
distributions.  For  example,  loading  variability  is  often  modeled  nonparametrically 
through  an  exceedance  curve  that  defines  the  number  of  occurrences  of  loads  (often 
accelerations  or  g-forces)  that  exceed  different  load  levels  [1,2].  The  upper  limit  to  the 
exceedance  curve  is  presumed,  but  often  not  verified,  to  be  sufficiently  large  such  that 
further  data  collection  will  not  affect  the  design  of  the  structure.  A  similar  condition 
arises  with  respect  to  defects  in  gas  turbine  materials  where  a  nonparametric  distribution 
is  used  to  model  the  probability  of  defect  size  [3].  This  distribution  is  then  used  as  part  of 
the  process  to  certify  that  the  design  is  within  allowable  limits  on  the  probability-of- 
failure.  A  truncated  distribution  in  the  context  of  this  research  encompasses  any 
distribution  that  has  non-zero  PDF  value  at  either  bound  of  its  domain. 
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Establishment  of  the  allowables  for  material  qualification  for  a  new  material  or  a 


new  application  of  the  material  is  typically  based  on  extensive  testing  and  is  often  costly 
in  terms  of  schedule  and  funds.  Therefore,  it  is  useful  to  have  some  tools  to  assess,  in  the 
design  stage,  the  impact  the  bounds  have  on  the  product  reliability,  cost, 
manufacturability,  etc.  Therefore,  methods  to  compute  the  sensitivity  of  the  reliability  (or 
probability-of-failure)  of  a  component  or  system  with  respect  to  the  bounds  of  truncated 
random  variables  would  be  a  useful  tool  to  assist  in  setting  appropriate  values  for  the 
bounds  and  to  assess  the  performance  of  the  design. 

The  finite  difference  method  can  be  used  to  estimate  the  numerical  derivative  of 
the  probability-of-failure  with  respect  to  a  random  variable  bound,  e.g.  perturb  a  bound, 

dPf  P(d  +  A9)-P(d) 


rerun  the  analysis  and  compute  the  ratio  of  the  differences 


36, 


A6 


- .  This 


approach  is  error  prone,  however,  as  the  amount  of  the  perturbation,  A6.  can  strongly 
affect  the  accuracy  of  the  calculation  and  the  numerical  estimate  usually  suffers  from 
subtraction  error.  Moreover,  a  small  perturbation  is  needed  to  obtain  an  accurate 
derivative.  As  a  result,  the  derivative  estimate  requires  a  subtraction  of  two  probabilities 
that  differ  only  slightly.  If  sampling  is  used  to  obtain  the  probability,  a  very  large  number 
of  samples  may  be  required  to  reduce  sampling  variance  before  the  subtraction  operation. 
These  issues  make  finite  differencing  error  prone  and  arduous  before  obtaining 
confidence  in  the  result. 

There  are  a  number  of  probabilistic  sensitivity  methods  that  provide  valuable 
information  during  the  design  process,  often  with  little  additional  computational  cost. 
Frey  and  S.R.  Path  provide  an  overview  article  discussing  ten  sensitivity  methods,  both 
probabilistic  and  deterministic  [4].  Helton  et  al.  discusses  the  use  of  scatter  plots, 
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regression  analysis,  and  Spearman  or  Pearson  correlation  and  other  methods  as 
qualitative  and  quantitative  metrics  for  sampling  methods  as  inexpensive  approaches  to 
discerning  the  contribution  of  the  variance  of  each  random  variable  to  the  output  [5]. 

Variance-based  methods  based  [6,  7]  are  capable  of  apportioning  the  amount  of 
variance  in  the  output  variance  to  the  variance  of  the  inputs,  thus,  providing  a  method  to 
rank  order  the  key  inputs.  The  method  is  applicable  to  monotonic  and  nonmonotonic 
models.  Calculation  of  the  sensitivity  indices  requires  computation  in  addition  to  the 
probability-of-failure  calculation  by  requiring  multiple  multidimensional  integrals.  The 
Fourier  Amplitude  Sensitivity  Test  (FAST)  method  can  be  used  to  reduce  the 
multidimensional  integrals  into  a  one-dimensional  integral  [8].  Lui  et  al.  use  the 
Kullback-Leibler  relative  entropy-based  method  to  evaluate  the  impact  of  a  random 
variable  on  a  design  performance  by  measuring  the  divergence  between  two  probability 
density  functions  of  the  performance  response,  obtained  before  and  after  the  variation 
reduction  of  the  random  variable  [9] .  The  application  is  similar  to  the  variance-based 
methods  but  is  not  limited  to  differences  in  the  second  moment. 

A  number  of  sensitivity  methods  are  available  for  the  First  Order  Reliability 
Method  (FORM).  Sensitivity  factors  (derivatives  of  the  safety  index  with  respect  to  the 
random  variables)  [10],  derivatives  of  the  probability  of  failure  with  respect  to  the 
random  variable  parameters,  e.g.,  dP  Id ji.  dP/do  [10]  and  omission  factors  [11]  are 
computed  as  by-products  of  an  analysis.  Each  of  these  indexes  provides  an  indication  as 
to  the  importance  of  the  parameters. 

Various  authors  develop  and  discuss  the  “Score  Function”  method  for  the 
computation  of  partial  derivatives  of  a  performance  function  (probability-of-failure  or 
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response  moment)  with  respect  to  parameters  of  the  underlying  input  probability 
distributions  [12-17].  This  method  provides  local  partial  derivatives  of  the  probability-of- 
failure  or  response  moments  with  respect  to  the  parameters  of  the  input  PDFs,  e.g., 
dP /dfl,  dP /do .  Implementation  of  the  methodology  is  convenient  using  sampling 
methods.  A  significant  advantage  is  that  negligible  additional  computing  time  is  required 
to  determine  the  sensitivities. 

None  of  aforementioned  methods  explicitly  treat  the  sensitivity  of  the  probability- 
of-failure  or  the  response  moments  with  respect  to  the  bounds  of  a  distribution. 

Therefore,  a  new  methodology  is  derived  and  demonstrated  below.  The  mathematical 
formulation  follows  the  development  of  the  Score  Function  method  [13].  However,  as 
shown  below,  an  additional  flux  integral  is  required  that  is  not  present  in  the  traditional 
Score  Function  formulation.  This  flux  integral  is  also  a  probability  integral.  On  the  other 
hand,  other  aspects  of  this  formulation  simplify  relative  to  the  traditional  Score  Function 
approach.  The  implications  are  that  the  resulting  equations  can  be  solved  with  any 
probabilistic  method  such  as  sampling,  the  First  Order  Reliability  Method,  conditional 
expectation,  importance  sampling,  and  others. 

The  rest  of  this  paper  is  organized  as  follows.  The  basic  formulation  of  the 
governing  equations  for  sensitivities  of  the  probability-of-failure  is  presented  in  Section 
2.  Section  3  describes  the  extension  of  the  method  to  compute  sensitivities  with  respect  to 
response  moments.  Section  4  provides  numerical  examples  and  Section  5  contains 
discussion  and  conclusions. 
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2.  Basic  Formulation 


2.1  Governing  Equations 

The  methodology  developed  here  can  be  considered  as  an  extension  of  the  Score 
Function  method  [13].  Given  a  joint  probability  density  function  (JPDF)  /X(x),  the 
probability-of-failure  can  be  defined  as 


P,  =  f  f  (x)dx 

f  J  g(x)z 0JxK 


(1) 


where  X  is  a  vector  of  random  variables  of  length  N,  X  =  (Xp. .  .,XN),  g  is  the  limit  state 
function  where  g(x)  <  0  defines  the  failure  domain,  and  the  integral  is  N  dimensional. 

In  this  paper,  the  focus  is  on  developing  sensitivities  with  respect  to  bounds  of 
truncated  distributions,  that  is,  probability  distributions  whose  range  is  finite  in  at  least 
one  direction  and  whose  PDF  value  at  the  finite  range  is  nonzero.  In  addition,  only 
rectangular  truncation  is  considered,  i.e.,  xi  =  6i :  6i  E  R"  where  0,.  corresponds  to  the 
lower  or  upper  bound  of  the  distribution  fx  .  Sensitivities  with  respect  to  means  and 

standard  deviations  of  truncated  distributions  have  been  discussed  previously  [18]. 
Related  developments  on  sensitivities  with  respect  to  PDF  parameters  for  nontruncated 
distributions  can  be  found  in  [13,18] 

Introducing  the  indicator  function  /(x),  which  is  defined  as  equal  to  one  if 
g(x)  <  0  and  zero  otherwise,  the  probability-of-failure  integral  can  be  written 


Pf  =  f  /(x)  fx(x)dX  (2) 
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The  derivative  of  the  probability  integral  with  respect  to  a  parameter  of  a  random 
variable  that  affects  the  boundary  can  be  determined  by  using  the  idea  of  the  classical 
Reynold’s  Transport  Theorem: 

—  f  W(x,t)dV  =  f^V(X,t)dV  +  f  W(x,t)  VjHjdS  (3) 

Dt  y  y  St  £ 

where  'P  is  a  property  of  the  continuum,  v  denotes  the  velocity  of  the  material,  n 
represents  the  unit  normal  along  the  boundary  S,  and  V  is  the  volume  enclosed  by  S. 
The  surface  integral  term  in  Eq.  (3)  is  the  value  of  XP  on  the  boundary  multiplied  by  the 
volume  swept  by  the  particles  on  the  boundary  in  the  time  interval  dt,  integrated  over 
dS.  This  term  can  be  considered  as  a  flux  of  the  property  'F  over  the  surface  S.  The  total 
derivative,  D/Dt,  is  also  known  as  the  material  derivative  [19]. 

The  concept  of  the  material  derivative  can  be  utilized  to  take  the  derivative  of  the 
probability  integral,  Eq.  (2),  with  respect  to  a  bound  of  a  random  variable  PDF.  Here,  the 
independent  parameter,  £f ,  is  a  bound  of  the  distribution  representing  Xf,  the  JPDF  is 
equivalent  to  the  property  XP,  the  volume  is  the  N  dimensional  random  variable  space, 
and  S  is  the  surface  of  the  random  variable  space  remaining  when  random  variable  Xj  is 
set  to  the  bound  6r  For  rectangular  truncation,  the  surface  S  is  straightforward  to 
compute  as  the  independent  parameter,  0  ,  is  a  bound  of  the  N  dimensional  random 
variable  space. 

The  unit  normal  and  the  equivalent  velocity  term  and  their  relation  can  be 
discerned  from  a  problem  of  two  random  variables,  see  Figure  1.  Since  the  independent 
parameter  d  is  an  element  of  X,  the  velocity  becomes  v  =  dxldd  =  1.  At  the  lower  bound 
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v  and  n  are  in  opposite  directions,  hence,  the  dot  product  v  -n  •  equals  -1.  At  the  upper 
bound  VjHj  equals  +1. 

Applying  these  concepts  to  the  probability-of-failure  integral  yields  the  equation 

HD  00  00 

=  f  /(x,flj)  KM*  ±  f  i(y,0,)fx  (y-d,)dy  (4) 

86 :  J  86,  J 


where  Y  denotes  the  vector  of  the  random  variables  X  but  excluding  Xn  i.e., 

Y  =  (A1,...,A;_1,A(+1,...,Aa,),  /x( y,0,.)  represents  the  conditional  JPDF  given  xj  =  dj,and 
the  “+”  sign  is  used  when  0i  represents  the  upper  bound  and  for  the  lower  bound. 


Using  the  concept  of  kernel  functions  [18]  defined  as  Kg  (X)  = 


#x(X)  1 

/x(X)’ 


Eq. 


(4)  can  be  written 


-jr  =  f  I(x)Ke(x)fx(x)dx±  f  /(y,0,)./x(y,0,)dy 

^  i  -oo  -oo 


(5) 


The  kernel  functions  with  respect  to  a  bound  for  a  truncation  value  of  random 
variable  A(  are  (see  appendix  A) 


K 


a 


#x(X)  1 

da  /x(X) 


=  fx,  («) 


(6) 


k:,. 


#x(x)  1 

db  /X(x) 


=  -/xlW 


(7) 


where  r/  and  b  represent  the  lower  and  upper  bound,  respectively,  of  A( .  It  is  quite 
remarkable  that  the  magnitude  of  the  kernel  functions  are  independent  of  X,  equal  to  the 
value  of  the  marginal  PDF  of  Xi  evaluated  at  the  respective  bound,  and  independent  of 
the  form  of  the  JPDF.  In  a  subsequent  section,  we  demonstrate  that  the  kernel  functions 
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must  have  these  mathematical  forms  in  order  to  satisfy  the  requirement  that  the  derivative 
of  the  probability-of-failure  equals  zero  for  the  case  where  the  indicator  function  is 
everywhere  one. 

oo 

The  flux  surface  integral  rf  =  J I(y,6)  fx(y,d)dy  can  be  rewritten  since  the  JPDF 

—oo 

can  be  simplified  as  /x(y,0)  =  ,fx  (0)/Y( y)  yielding 


if  =  fx,  (0)  f  I(y,d)fY(y)dy  =  fXt  (d)P(;  (8) 


where  Pf  =  f  I(y,6)fY(y)dy  and  the  subscript  /  denotes  the  dependence  on  the  indicator 

—oo 

function.  Since  Pf  is  a  probability  integral,  it  is  always  positive.  Using  the  kernel 
functions  from  Eqs.  (6)  and  (7),  the  surface  integral,  Eq.  (8),  the  fact  that  Pf  =  £(/(x)|, 
and  the  proper  sign  for  each  bound,  the  equations  for  the  sensitivities  become 


c7Pl 

da. 


fx  ( a)(pf-p “) 


(9) 


dPf 

db. 


■  ~fx  (b)(P f  ~  P,) 


(10) 


Thus,  the  effort  to  obtain  the  sensitivities  becomes  one  of  evaluating  the  flux  integrals  P° 
and  Pj  in  addition  to  the  probability-of-failure. 


2.2  Evaluation  of  flux  integral 

The  flux  integral  can  be  evaluated  by  integrating  the  JPDF  times  the  indicator 
function  over  the  surface  S  defined  by  the  condition  xt  -  0.  That  is,  one  integrates  the 
JPDF  in  the  failure  region  over  the  surface  S.  Thus,  the  dimension  of  the  flux  integral  is 
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N-l.  A  schematic  of  the  flux  in  two  dimensions  (two  random  variable  problem)  is  shown 
in  Figure  2  where  the  value  of  the  joint  PDF  along  the  upper  bound  of  Xt  is  outlined. 
Mathematically,  the  flux  is  computed  as 

r»=fx(e)po  (ID 

For  problems  with  one  random  variable,  the  flux  reduces  to  fx  {a)  over  the  lower  bound 
and  /x  {b)  over  the  upper  bound  and  Eqs.  (9)  and  (10)  reduce  to  the  well-known  Leibnitz 
rule  in  one  dimension. 

Figure  3  shows  the  failure  domain  for  several  problems  of  two  random  variables. 
In  each  case,  the  failure  domain  fully  encompasses  the  upper  bound  of  random  variable 
Xj  and  does  not  encompass  any  of  the  lower  bound.  Therefore,  given  the  same  JPDF, 
each  problem  has  the  same  flux  integral  with  respect  to  the  bounds,  that  is,  the  flux 
integral  depends  only  upon  the  failure  definition  along  the  boundary  S  and  not  by  the 
definition  of  failure  throughout  the  failure  domain.  For  the  examples  given  in  Figure  3, 

Pc;  =0  and  P,h  =  1 . 

Evaluation  of  Eq.  (11)  requires  an  N-l  dimensional  integral.  However,  the  flux 
integral  is  also  a  probability  integral,  which  can  be  considered  as  a  special  case  of  the 
probability-of-failure  integral,  Eq.  (2).  Therefore,  the  flux  integral  can  be  computed  using 
standard  probability  methods  such  as  Monte  Carlo  sampling,  First  Order  Reliability 
Method,  importance  sampling,  and  different  methods  can  be  used  to  compute  the 
probability-of-failure  and  the  fluxes. 

The  flux  integral  is  evaluated  by  integrating  the  JPDF  within  the  failure  domain 
over  the  surface  defined  by  xi  =  Ui.  However,  the  integration  may  proceed  over  the  safe 
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domain  defined  by  Pj  =  J  I(y,d)fY(y)dy,  where  /( y,d)  equals  one  in  the  safe  domain, 

—oo 

zero  otherwise  (the  complement  to  7(y,0)).  Then  the  flux  integral  can  be  computed 
indirectly  since  Pf  +  Pj  =1.  Then  Pf  =1  -  Pj  ■ 

Comparison  between  Eqs.  (11)  and  (2)  indicates  that  the  numerical  values  of  the 
probability-of-failure  and  the  flux  may  be  significantly  different  as  the  integrals  are  of 
different  dimensions  and  domains.  For  example,  consider  upper  right  graphic  in  Figure  3 
where  the  failure  domain  is  a  small  sliver  near  the  upper  bound  of  Xl .  As  a  result,  the 
probability-of-failure,  Pf ,  will  be  very  small  whereas  the  flux  integral  over  the  surface, 
Pjb,  will  be  equal  to  one  since  the  integration  of  the  joint  PDF  with  respect  to  the 
remaining  random  variables  integrates  to  one.  Thus,  Pf  and  Pb  will  have  significantly 
different  values. 

The  flux  integrals  are  particularly  straightforward  to  compute  when  using 
sampling  as  they  are  a  subset  of  the  probability  integral.  The  sampling  code  for  the 
probability-of-failure  can  be  used  with  the  only  change  of  generating  samples  at  the 
appropriate  bound  for  the  affected  random  variable.  The  necessary  equations  are 

I  m 

p, <12) 

I  m 

r?—Z«j„e>  (13) 

mc! 

If  sampling  is  used  to  compute  Pf ,  it  is  possible  to  reuse  the  samples  to  compute 

P°  by  projecting  the  samples  onto  the  surface  Xt  =  a  or  X t  =  b  thereby  obtaining  the 
sensitivities  at  negligible  computational  cost. 
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The  flux  integrals  can  also  be  computed  using  First/Second  Order  Reliability 
Method  (FORM/SORM)  by  recasting  the  integrals  in  the  N-lth  dimension:  the  limit  state 
is  modified  to  g(x  I  xt  =  0)  =  0,  variable  Xj  is  no  longer  random,  and  fx  =  fx  (0). 


2.3  Maximum  and  Minimum  Sensitivities 


The  probabilistic  sensitivities,  defined  in  Eqs.  (9)  and  (10),  are  dependent  upon 
the  probability  integrals  Pf  and  P/.  Given  that  probability  integrals  are  bounded  by  zero 

and  one,  it  is  straightforward  to  estimate  the  maximum  and  minimum  magnitude  the 
sensitivities  may  obtain.  These  estimates  may  be  sufficient  for  preliminary  design 
purposes. 

The  largest  positive  values  of  the  sensitivities  occur  when  P‘‘  =0,  and  P/  =  1, 
yielding 


dPf 


dat 
dP  t 


db; 


=  fXi(a)Pf 
=  fx(b)(\-Pf) 


(14) 


Conversely,  the  minimum  values  of  the  sensitivities  occur  when  P“  =1,  and  P/’  =0, 
yielding 


dP, 


dai 

dPr 


=  -fx(a)(l-Pf) 


db: 


Neg 


Neg 


- fx.(b)Pf 


(15) 


In  certain  instances,  it  may  be  possible  to  have  a  priori  knowledge  of  the  failure 
domain  along  a  bound  of  the  truncated  random  variable,  and,  therefore,  have  an  estimate 
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of  the  flux  integral  and  the  sensitivity.  For  example,  if  the  random  variable  represents  the 
load  on  a  component,  it  may  be  known  that  the  component  will  fail  for  a  load  equal  to  the 
upper  bound  of  the  random  variable  and  not  fail  for  a  load  equal  to  the  lower  bound  of  the 
random  variable,  regardless  of  the  values  of  the  other  random  variables.  In  such  a  case, 
the  failure  domain  along  the  surfaces  xt  =  a  and  xi  -  b  is  known  a  priori  and  in  this  case 
encompasses  the  entire  surface  defined  by  xt  =  b  and  none  of  the  surface  defined  by 
Xj  -  a.  Along  the  lower  bound  the  flux  integral  is  zero  and  along  the  upper  bound  the 
flux  integral  is  one.  In  this  case,  the  sensitivities  are  given  in  Eqs.  (14)  and  no  explicit 
evaluation  of  the  flux  integrals  is  required. 

2.4  Properties  of  the  Sensitivities 

It  is  known  that  the  sensitivities  must  satisfy  certain  properties  [13,  18].  In 
particular,  if  the  failure  domain  encompasses  the  entire  random  variable  domain  7(x)  =  1 
everywhere,  then  the  probability-of-failure  is  always  one  regardless  of  the  values  of  the 
parameters;  therefore,  dPf  IdO  =  0.  This  condition  must  be  satisfied  with  respect  to  any 
parameter  of  the  PDFs  including  bounds.  When  the  failure  domain  encompasses  the 
entire  random  variable  domain  Pf  =  P°  =  1  and  the  sensitivities  are  dPf  /dd  =  0.  From 
these  results,  we  see  that  the  kernel  functions  must  have  the  forms  Ka  =  fx  {a), 

Ka  =  -fx  ( a )  regardless  of  distribution  type  in  order  for  dPf  IdQ-  0  to  be  satisfied. 
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2.5  Variance  estimates 


Variance  estimates  for  dPf  /dd  can  be  obtained  by  computing  the  variance  of  Eqs. 
(9)  and  (10).  The  formulation  is  straightforward  as  the  variance  of  a  linear  function  of 
random  variables  X,  e.g.,  Y  =  '^AjXj,  where  A  are  constants,  is  known  analytically  as 
[20] 

Var(Y)  =  22AiAJPuoiaj  (16) 

>=i  j= i 

where  cr  represents  the  standard  deviation  of  Xn  and  pi}  denotes  the  correlation  between 
Xt  and  X j.  Applying  Eq.  (16)  to  Eqs.  (9)  and  (10),  yields 

dP 

V\r^\  =  fl  (6)(V[Pf]  +  V[P?]  -  2  po[PfMP?i)  (17) 

where  p  is  the  correlation  between  Pf  and  P,0.  It  is  clear  from  Eq.  (17)  that  the 
correlation  between  Pf  and  P0  affects  the  variance  of  the  sensitivity  estimate. 
Independent  sampling  can  be  chosen  to  ensure  p  «  0.  A  better  approach  is  to  employ 
“common”  variables  such  that  the  correlation  is  greater  than  zero,  i.e.,  p  >  0,  to  minimize 
the  variance  [21,  20].  This  is  easily  accomplished  by  using  the  same  samples  for  the  Pf 
and  Pj  integrals.  The  reuse  of  samples  by  projecting  the  samples  onto  the  bound  of 
interest  ensures  positive  correlation  and  provides  an  estimate  of  the  flux  integral  without 
additional  computational  cost. 
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The  variance  estimates  of  Pf  and  P1/  are  well  known  if  Monte  Carlo  sampling  is 


N 


used  to  compute  the  probabilities,  e.g.,  V[Pf] 

i  \ 

Pf  and  P,a  are  the  sampling  estimates  of  P ,  and  P“ . 


and  V[Pf\ 


p;‘(\  -  p;‘) 
N 


where 


3.  Extension  to  Response  Moments 

In  many  mathematical  and  engineering  problems,  the  sensitivity  of  the  response 
moments  (mean,  /j_ ,  and  standard  deviation,  a.)  are  of  interest,  where  Z(X)  is  an 
arbitrary  output  function  of  the  random  variables.  The  methodology  derived  in  the 
previous  sections  can  be  extended  to  provide  the  probabilistic  sensitivities  of  the  response 
moments  with  respect  to  the  bounds,  e.g.,  dflz  Ida  or  doz  Ida ,  see  Appendix  B. 


3.1  Probabilistic  sensitivity  of  the  mean  response  with  respect  to 
bounds 


The  sensitivities  of  the  response  mean,  jiz , 
(see  Appendix  B) 

d-i±  =  fx\a)(nz-PZ) 

dat 


db, 


-fX;  (a)(pz  ~  Pz  ) 


to  the  bounds  of  a  truncated  PDF  are, 


(18) 


The  flux  calculation  here  is  modified  from  Eq.  (1 1)  in  that  the  indicator  function 
is  now  replaced  by  the  response  function;  therefore,  the  flux  can  be  written  as 


T|  =  /  Z(y,d)fx(y,6)dy  =  fx_  ( 6)Pez  (19) 
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where 


Pez  =  f  Z(y,d)fY(y,d)dy  (20) 

—oo 

The  integration  domain  is  again  defined  by  the  surface  S  defined  by  xt  =  Uj. 
Although  the  probability  integral  must  be  positive,  the  flux  integral  P®  can  be  either 
positive  or  negative  depending  upon  the  response  function  Z. 

3.2  Probabilistic  sensitivity  of  the  standard  deviation  of  the  response 
with  respect  to  bounds 

The  sensitivities  of  the  response  standard  deviation,  oz  to  the  bounds  of  a 
truncated  PDF  are,  (see  Appendix  B) 

^T  =  fxXa){(?z  -A4  +2pzPl-Paz2}l{2oz)  (21) 

da  '  z 

=  -fXi  m<4  -  +  2^zfz"  -  f”. )  /(2oz)  (22) 

do  '  z 

oo 

where  P^  =  J Z2(y,d)fY(y)dy.  P°zl  must  be  always  positive. 

—oo 

4.  Numerical  Examples 

A  two  random  variable  problem  amenable  to  exact  integration  both  component 
and  system,  is  solved  in  detail  to  demonstrate  the  methodology.  This  problem  also 
provides  a  means  to  generate  plots  that  clearly  show  the  flux  integrals. 
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4.1  Two  random  variable  problem 

A  two  dimensional  example  is  solved  to  illustrate  the  concepts.  Consider  a  limit 
state  of  g(r,s )  =  r-  s  with  R  and  S  modeled  as  independent  random  variables.  The 
indicator  function  defines  the  failure  region  as 

I[r,s]  =  1  if  g(r,s)  <  0 

0  otherwise 

R  is  a  standard  normal  distribution  and  S  is  a  uniform  distribution  with  bounds  a  =  0  and 
b  =  1. 


/*(r)  = 


Exp[-r2  /2] 


—  oo  <  r  <  oo 


(24) 


fs(s)  =  l  0<5<1 
0  otherwise 


(25) 


4.2  Probability-of-failure  sensitivities 

The  JPDF  over  a  portion  of  the  failure  domain  ( -4  <  r  <  4;  0  <  s  <  1)  is  shown  in 
Figure  4  with  a  direct  view  of  the  face  5  =  1.  The  exact  solution  from  integration  is 
Pf  =  0.68437.  The  JPDF  with  a  direct  view  of  the  face  S  =  0  is  shown  in  Figure  5.  The 
fluxes  can  be  evaluated  using  one-dimensional  integrals 

oo 

P,a=fs(a)jl(r,a)fR(r)dr  (26) 


Pt=fs(b)$I(r,b)fR(r)dr  (27) 
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Using  symbolic  integration,  P“  =0.5  and  P]'  =0.8413.  From  Figures  4  and  5,  it  is  clear 
that  the  integral  of  fR  over  the  failure  surface  where  5=0  equals  FR  (0)  =  <f>(0)  and  the 
integral  of  f R  over  the  failure  surface  where  5=1  equals  PR(\)  =  $(1) . 

Table  1  summarizes  results  obtained  using  both  the  flux-based  methodology  and 
the  standard  finite  difference  (forward  differencing)  method.  The  derivatives  of  the 
probability-of-failure  with  respect  to  the  bounds  were  computed  using  Eqs.  (9)  and  (10). 
Finite  difference  estimates  were  computing  using  symbolic  integration  with  a  forward 
step  size  of  0.00001.  The  minimum  and  maximum  sensitivities  are 
dPf  Ida  =  [-0.3157,0.6843],  and  dPf  tdb  =  [-0.6843,0.3157] . 


1  N 

Sampling  estimates  were  computed  using  the  equations  P.  =  —  V  I(rk,sk), 

^  J  A I 


1  N  j  N 

P‘‘  =  — ^  I(rk,a),  and  Pk  ~  — ^  /(/;,/?)  with  common  random  variables  [21].  FORM 


analysis  can  also  be  used  to  compute  the  flux  integrals.  In  this  example,  FORM 
reproduces  the  exact  results. 

The  sensitivities  for  the  flux-based  and  finite  difference  methods  using  exact 
integration  are  very  close,  as  expected.  The  results  using  sampling  are  also  in  good 
agreement;  however,  note  that  104  samples  were  used  for  the  flux-based  approach  to 
obtain  a  solution  with  good  accuracy  versus  106  samples  for  the  finite  difference 
sampling-based  approach.  The  superiority  of  the  flux-based  approach  compared  to  the 
standard  finite  difference  (forward  differencing)  approach  using  sampling  can  be  shown 
clearly  by  comparing  the  95%  confidence  bounds  and  the  coefficient  of  variation  (COV  = 
standard  deviation/mean)  obtained  using  both  methods  obtained  from  100  trials.  Figures 
6a  and  6b  show  plots  of  the  95%  confidence  bounds  for  finite  difference  (dashed)  versus 
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flux-based  (solid)  for  dPf/da  and  dPf  /db,  respectively.  The  bounds  for  the  flux-based 

approach  are  so  much  narrower  than  the  finite  difference  approach  that  they  show  almost 
as  a  straight  line.  The  bounds  are  so  wide  for  the  finite  difference  method  that  any 
solution  obtained  is  completely  unreliable  until  the  number  of  samples  approaches  one 
million. 

A  closer  examination  of  the  benefits  of  using  common  random  variables  during 
the  sampling  process  is  provided  in  Table  2.  Variance  results  for  negative,  approximately 
zero,  and  positive  sampling  correlation  are  provided.  The  results  clearly  show  that 
positive  correlation  provides  approximately  a  three  times  reduction  in  the  variance  of  the 
sensitivities  with  respect  to  independent  sampling;  at  zero  cost.  Positive  correlation  was 
accomplished  simply  by  using  the  same  samples  to  compute  the  probability-of-failure 
and  the  fluxes.  All  subsequent  results  using  sampling  shown  in  the  paper  were  computed 
using  common  variables,  thereby  inducing  positive  correlation. 

Table  3  shows  the  coefficient  of  variation  of  the  two  methods  as  a  function  of  the 
number  of  samples.  The  COV  for  the  finite  difference  method  is  approximately  two 
orders  of  magnitude  larger  than  that  obtained  using  the  flux-based  method  for  the  same 
number  of  samples.  These  results  imply  that  there  are  approximately  4  orders  of 
magnitude  difference  in  the  number  of  samples  required  to  achieve  similar  accuracy 
using  both  methods.  The  explanation  is  that  the  finite  difference  method  requires  an 
approximation  of  a  limiting  process  estimated  by  subtracting  two  near-equal  numbers. 
The  flux-based  approach,  on  the  other  hand,  requires  no  limiting  process  nor  subtraction 
of  near-equal  numbers.  The  results  from  Table  3  indicate  that  using  the  flux-based 
approach  with  103  samples  is  superior  to  the  finite  difference  method  with  106  samples. 


21 


4.3  Response  moments 


Consider  the  previous  PDFs  with  a  response  function 
Z(r,s)  =  5  +  5r-2r-5s-As2.  Plots  of  Z  times  the  JPDF  are  shown  in  Figures  7a  and  7b 
and  Z2  times  the  JPDF  are  shown  in  Figures  8a  and  8b.  The  profiles  of  the  flux  of  Z  fR 
and  Z2  fR  along  the  boundaries  s  =  a  and  s  =  b  are  clearly  evident.  Using  numerical 
integration,  the  moments  of  Z  are  jlz  =-0.8333  and  oz  =6.311.  The  fluxes  are 
rz(y,a)  =  3,  rz(y,Z?)  =  -6,  rz2(y,a)  =  42,  and  rz,(y,^)  =  69.  The  sensitivities  are  given 

in  Table  4.  The  flux-based  and  finite  difference  results  using  integration  were  very  close, 
as  expected. 

The  expected  values  of  the  sensitivities  using  sampling  showed  close  agreement 
with  the  numerical  integration  results.  However,  similar  to  the  sensitivities  with  respect  to 
probability-of-failure  shown  in  Tables  1  through  3,  the  COV  of  the  sensitivities  of  the 
response  moments  using  the  finite  difference  approach  was  much  larger  than  the  flux- 
based  method;  10%  for  djlz  Ida  using  finite  difference  with  106  samples  versus  0.6%  for 
the  flux-based  approach  using  104  samples.  The  COV  results  for  doz  Ida  were  50%  with 
for  finite  difference  with  106  samples  and  5%  for  the  flux-based  approach  with  104 
samples.  Similar  results  were  seen  for  djlz  I db  and  doz  I db. 

4.4.  System  reliability 

The  previous  example  was  modified  to  demonstrate  the  applicability  of  the 
sensitivity  methodology  to  system  reliability  problems.  Two  limit  states  were  considered 
g{(r,s)  =  r-  s  and  g2(r,s)  =  5  -5r  +  2r2  +  5s  +  As2.  The  system  was  analyzed  as  both  a 
series  and  a  parallel  system. 
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4.4.1  Series  system 


For  the  series  system,  failure  was  defined  as  Pf  =  P[g,  <  0  U  <?2  -  0  ]  •  The  JPDF 

in  the  failure  region  is  shown  in  Figure  9.  The  JPDF  along  the  bounds  a  and  b  needed  for 
the  flux  calculations  are  shown  in  Figure  10a  and  Figure  10b,  respectively.  The  results 
were  Pf  =0.70375  (exact  integration),  P‘‘  =  <I>(0)  +  «T>(— 3. 26557)  =  0.5000546,  Pb  =1, 
dPf  Ida  =  0.203204  and  dPf  Idb  =  0.29625 .  The  sensitivities  are  shown  in  Table  5.  The 
minimum  and  maximum  sensitivities  are  dPf  Ida  =  [-0.2962,0.7038],  and 
dPf  Idb  =  [-0.7038,0.2963] .  As  shown  in  previous  examples,  the  COV  estimates  using 

sampling  for  the  finite  difference  method  with  106  samples  was  approximately  an  order 
of  magnitude  larger  than  the  COV  for  the  flux-based  method  using  104  samples. 

4.4.2  Parallel  system 

For  the  parallel  system,  failure  was  defined  as  Pf  =  P[g,  <  0  P|  g2  <  0  ] .  The  JPDF 
in  the  failure  region  is  shown  in  Figure  11.  The  JPDF  along  the  bound  S  =  0  is  shown  in 
Figure  12  and  it  is  clear  from  the  Figure  1 1  that  the  flux  across  the  bound  S  =  1  is  zero. 
The  results  from  numerical  integration  were  Pf  =  0.233742, 

P;  =<F(0)-  $(-0.76556)  =  0.278031,  P,b  =0,  dPf  Ida  =  -0.0442901  and 

dPf  Idb  =  -0.233742 .  The  minimum  and  maximum  sensitivities  are 

dPf  Ida  =  [-0.7663,0.2337],  and  dPf  Idb  =  [-0.2337,0.7663] .  The  sensitivities  are  shown 

in  Table  6.  Again,  the  COV  for  the  flux-based  method  was  significantly  smaller  than 
finite  difference  even  with  fewer  samples. 
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5.  Discussion  and  Conclusions 


Efficient  evaluation  of  the  sensitivity  of  the  probability-of-failure  or  the  response 
moments  to  the  bounds  of  truncated  distributions  can  provide  useful  information  in  the 
design  stage  in  order  to  optimize  product  reliability,  minimize  cost,  determine  quality 
assurance  procedures,  etc.  The  method  outlined  here  can  be  used  to  compute  these 
sensitivities  with  a  significant  improvement  in  computational  efficiency  over  standard 
finite  difference  methods. 

The  methodology  is  based  upon  an  application  of  the  material  derivative  concept 
to  the  probability-of-failure  or  the  response  moment  integrals  thus  yielding  a  flux  integral 
that  must  be  computed  in  addition  to  the  standard  probability-of-failure  or  response 
integrals.  The  methodology  is  applicable  to  any  limit  state  formulation,  either  component 
or  system,  and  any  random  variables  described  by  a  truncated  joint  probability  density 
function  containing  either  correlated,  e.g.,  truncated  multivariate  normal,  or  independent 
random  variables. 

The  sensitivities  require  a  supplemental  flux  integral  for  each  bound  that,  when 
combined  with  the  probability-of-failure  and  kernel  functions,  provides  the  sensitivity 
with  respect  to  the  bound  of  a  truncated  distribution.  However,  the  flux  integral  itself  is  a 
probability  integral  and,  therefore,  amenable  to  solution  using  existing  probabilistic 
methods.  Simple  approximations  to  the  flux  integral  are  easily  discerned  which  may 
avoid  explicit  calculation  of  the  integral  and  provide  an  upper  bound  on  the  sensitivities. 
The  magnitudes  of  the  probability-of-failure  and  the  flux  integrals  maybe  significantly 
different  and  different  methods  may  be  used  to  compute  each. 


24 


If  sampling  is  used  to  compute  the  probability-of-failure,  the  samples  can  be 
reused  to  compute  flux  integrals  by  projecting  the  samples  to  the  bound  of  interest.  Thus, 
the  flux  integrals  are  computed  with  negligible  additional  computation. 

The  superiority  of  the  flux-based  approach  over  the  standard  finite  difference 
method  was  clearly  evident  from  numerical  studies  using  Monte  Carlo  sampling  that 
indicated  that  the  estimate  of  the  sensitivities  using  the  flux-integral  approach  required 
approximately  4  orders  of  magnitude  fewer  samples  for  the  same  accuracy  as  a  standard 
finite  difference  approach. 
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Appendices 


Appendix  A  -  Kernel  functions  with  respect  to  PDF  bounds 

The  JPDF  subjected  to  rectangular  truncation,  i.e.,  xt  =  dl :  6i  E  R"  where  0 
corresponds  to  the  lower  (a)  or  upper  ( b )  bound,  can  be  written 


/x(X) 


/x(X) 

Fx(b)-Fx(a) 


0 


a  <  x,  <  b 
otherwise 


(A-l) 


where  /X(x)  is  the  non-truncated  JPDF,  a  and  b  are  the  lower  and  upper  bounds  of 
random  variable  Xn  and  Fx  denotes  the  CDF  of  Xr 
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Consider  the  derivative  with  respect  to  the  upper  bound  b 


ffx(x)  _  d  I  /X(x) 
db  db\FXi(b)-FXi(a)) 


(  „  ^ 

-/x(X) 

dFx(b) 

(FSi{b)-FSi(ajf  t 

db 

{  -A(X) 

fxib)  = 

(N 

r  F 

1 

ST 

-/x(X) 

fXi 

(A-2) 


Fx(b)-FXi(a)Fx(b)-FXi(a) 

-fx(X)fx(b) 


Therefore,  the  kernel  function  with  respect  to  a  rectangular  bound  xt  =  b,  where  b  is  an 
upper  bound,  becomes 


j/x(X)  1  = 

db  /x(X)  •'*' 


(A-3) 


Similarly 


K 


a 


tfxW  1 
da  /x(X) 


=  /x,  («) 


(A-4) 


That  is,  the  magnitude  of  the  kernel  function  with  respect  to  a  bound  of  a  random  variable 
X t  is  equal  to  the  marginal  distribution  of  X t  evaluated  at  the  bound.  The  sign  of  the 
kernel  function  is  dependent  on  the  bound  in  consideration.  These  results  are  independent 
of  distribution  type. 
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Appendix  B  -  Derivation  of  moment  sensitivities 


The  sensitivity  of  the  response  mean  jlz  with  respect  to  the  distribution  bounds  is 
derived  as  an  extension  of  the  approach  by  Wu  and  Mohanty  2005  by  adding  the  flux 
integral, 


oo  oo 

jt  =  ~  f  Z(x)fx(x)dx±  rz(y,0)  =  f  Z(x)Kg(x)fx(x)dx±  Tz(y,0)  = 

U  U-  U  V  j  _oo  _oo 

E{Z(x)Ke  (x)]  ±  Tz  (y,  0) 


where  Z(x)  represents  the  deterministic  model  response,  Ke  (x)  represents  the  kernel 
functions.  Using  the  kernel  functions  defined  in  Eqs.  (6)  and  (7), 


d^  =  fx(a)ilz-rz(y,a) 
da 

%  =  -A.(i>)fe  +  rz(y,4) 

db 


(B-2) 


Since  Tz(y,d)  =  fx{9)Pz,  the  sensitivities  become 


^  =  /XiW(fe-Pz“) 

da 


(B-3) 


The  sensitivities  of  the  standard  deviation  of  the  response  with  respect  to  the 
bounds  of  a  random  variable  PDF  can  be  determined  by  extending  the  derivation  of  Sues 
and  Cesare  [17]  through  the  addition  of  the  flux  integrals  to  yield 
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~rL~  =  (E\Z2Ka\ -  r2  (y  ,a)  -  2juz(E[ZkJ  -  Tz(y,a)))/(2oz)  = 
da  z 

(fx  Xa)E[Z2] -rz2(y,a) -2nz(fXi(a)E[Z] -Tz(y,a)))/(2oz)  = 
(fxXa)(crz+^)-fxXa)P^  -2nz(fXi(a)ilz  -  fxXa)P“))/(2oz)  = 
fxXa){ ce  -  jly  + 2 jizPz  -  P“, } / (2oz ) 

Where  we  have  used  the  fact  that  E[Z2]  =  crz+  jlz,  Tz(y,0)  =  fx{d)Pz  ,  and 
rz2(y,e)  =  fx(d)Pez2. 

A  similar  derivation  for  doz/db  yields 

^  =  (E[Z2Kh]  +  Tz2(y,b)-2^z{E[ZKh]  +  Tz(y,b)))X2az)  = 

(lb  z 

(-/* Xb)E[Z2]  +  fx Xb)Pbz2  - 2nz(-fxXb)E[Z\  +  fxXb)P*))/(2oz )  = 

+£)  +  fxXb)P‘,  -2/1 z(-fx(b)Hz  +  f\  (b)Py  ))b2oy }  = 

-A,  H°i  - #4  +  -  E-  M2°z) 
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Tables 


Flux-based 

Finite  Difference 

Integration 

Sampling 

Integration 

Sampling 

(exact) 

(104  samples) 

(exact) 

(106  samples) 

dpf 

0.1844 

0.18421 

0.1844 

.18761 

da 

dpf 

0.1560 

0. 157 1 1 

0.1570 

.16081 

db 

'expected  value  of  100  trials) 

Table  1  Probabilistic  sensitivity  results  for  limit  state  g  =  r  -  5  (R  standard  normal,  S 
uniform) 
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#  Samples 

3pi 

dpi 

dpf 

dpf 

dpf 

dpf 

da 

da 

da 

3b 

3b 

3b 

(P<  0) 

(p-  0) 

(P>  0) 

(P<  0) 

(P-  0) 

(P>  0) 

103 

8.81E-4 

4.64E-4 

1.37E-4 

4.84E-4 

3.48E-4 

1.29E-4 

(p=-0.75) 

(P~0) 

(p=0.82) 

(p=-0.46) 

(p~o) 

(p=0.78) 

104 

7.37E-4 

4.01E-5 

1.44E-5 

4.03E-4 

2.70E-5 

1.22E-5 

(p=-0.66) 

(p=0.09) 

(p=0.68) 

(p=-0.27) 

(p=0.09) 

(p=0.64) 

Table  2  Variance  of  sensitivity  estimates  dPf  1 36  as  a  function  of  sampling 
correlation  (100  trials)  -  actual  correlation  in  parentheses 
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Flux-based 

Finite  Difference 

Integration 

Sampling 

Integration 

Sampling 

(numerical) 

(104  samples) 

(numerical) 

(106  samples) 

-3.833 

-3.8391 

-3.833 

-3.8481 

da 

(0.6%) 

(10%) 

-5.167 

-5.161 1 

-5.167 

-5.1531 

db 

(0.5%) 

(8%) 

doz 

-0.6223 

-0.62801 

-0.6223 

-0.66911 

da 

(5%) 

(50%) 

doz 

1.573 

1.5691 

1.573 

1 .55 1 1 

db 

(2.5%) 

(22%) 

'expected  value  of  100  trials 


Table  4  Sensitivity  results  for  response  moments;  COV  results  in  parentheses  based 
on  100  trials 
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Flux-based 

Finite  Difference 

Integration 

Sampling 

Integration 

Sampling 

(exact) 

(104  samples) 

(exact) 

(106  samples) 

dpf 

0.2032 

0.203 11 

0.2032 

0.21161 

da 

(1.6%) 

(18%) 

dPf 

0.2963 

0.29631 

0.2965 

0.30161 

db 

(1.3%) 

(13%) 

'expected  value  of  100  trials 

Table  5  Sensitivity  results  for  series  system;  COV  results  in  parentheses  based  on 
100  trials 
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Flux-based 

Finite  Difference 

Integration 

Sampling 

Integration 

Sampling 

(exact) 

(104  samples) 

(exact) 

(106  samples) 

dpf 

-0.04429 

-.043961 

-0.04429 

-0.04821 1 

da 

(12%) 

(111%) 

dPf 

-0.2337 

-0.23311 

-0.2337 

-0.22071 

db 

(2%) 

(22%) 

expected  value  of  100  trials 

Table  6  Sensitivity  results  for  parallel  system;  COV  results  in  parentheses  based  on 
100  trials 
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Figure  2  Flux  of  the  JPDF  in  the  failure  region  over  the  surface  S 
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Figure  3  Two  dimensional  problems  with  the  same  flux  values  with  respect  to  lower 
and  upper  bounds 
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Figure  4  Joint  PDF  over  failure  domain  (direct  view  of  the  face  of  S  =  1) 
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Figure  5  Joint  PDF  over  failure  domain  (direct  view  of  the  face  of  5=0) 
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Figure  6a  95%  confidence  limits  (100  trials)  for  dPf  Ida  as  a  function  of  the  number 
of  samples  (dashed  -  finite  difference;  solid  -  flux  based) 
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Figure  6b  95%  confidence  limits  (100  trials)  for  dPf  idb  as  a  function  of  the  number 
of  samples  (dashed  -  finite  difference;  solid  -  flux  based) 
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Figure  7a  Plot  of  JPDF  times  Z  (direct  view  of  the  face  5=1) 
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1.0 


Figure  7b  Plot  of  JPDF  times  Z  (direct  view  of  the  face  5=0) 
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Figure  8a  Plot  of  JPDF  times  Z2  (direct  view  of  the  face  5=0) 
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Figure  8b  Plot  of  JPDF  times  Z2  (direct  view  of  the  face  S  =  1) 
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Figure  9  JPDF  of  series  system 
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Figure  10a  Flux  along  bound  S  =  0  for  series  system 
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Figure  11  JPDF  of  parallel  system 
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Figure  12  Flux  along  bound  a  of  JPDF  along  bound  a  for  parallel  system 
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